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. ABSTRACT 

' A Quadratic Maximum Likelihood (QML) estimator is applied to the WMAP 5 year 

low resolution maps to compute the CMB angular power spectra at large scales for 
both temperature and polarization. Estimates and error bars for the six angular power 

(N ■ 

spectra are provided up to € = 32 and compared, when possible, to those obtained by 
the WMAP team, without finding any inconsistency. The conditional likelihood slices 
are also computed for the Ci of all the six power spectra from £ = 2 to 10 through 
a pixel based likelihood code. Both the codes treat the covariance for (T, Q,U) in 
. a single matrix without employing any approximation. The inputs of both the codes 

(foreground reduced maps, related covariances and masks) are provided by the WMAP 
team. The peaks of the likelihood slices are always consistent with the QML estimates 
within the error bars, however an excellent agreement occurs when the QML estimates 
are used as fiducial power spectrum instead of the best-fit theoretical power spectrum. 
By the full computation of the conditional likelihood on the estimated spectra, the 
value of the temperature quadrupole is found to be less than 2 a away from the 
WMAP 5 yrs ACDM best-fit value. The BB spectrum is found well consistent with 
zero and upper limits on the B-modes are provided. The parity odd signals TB and 
' EB are found consistent with zero. 

Key words: cosmic microwave background - cosmology: theory - methods: numerical 
- methods: statistical - cosmology: observations 
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1 INTRODUCTION 

The Cosmic Microwave Background (CMB hereafter) is a 
powerful tool for investigating the properties of the early 
and present Universe. Under hypothesis of Gaussianity 
and statistical isotropy, the main cosmological information 
coming from a CMB map is contained in its temperature 
and polarization angular power spectra (APS hereafter). 
In recent years several experiments, with different de- 
tection techniques, have measured the CMB anis o tropie s 
with increas i ngly b e tter precision (iHinshaw et al. ll2003h: 
iKogut et ai] l|2003l ): Ijones et al.l (|2006l ): iPiacentini et al] 
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Many methods have been developed to give un- 
bias ed estimates of the CMB power spectra a s Mas - 
ter jHivon et alll2002l), Cross-Spe c tra (ISaha et al.l l|2006l ). 
iPolenta et all ll2005h. ICrain et all l|2009l )). At high multi- 
poles {e > 30, 'Efst athioul (|2004| )) the so called pseudo - Ci 
algorithms are usually preferred to others techniques. These 
methods in fact implement the estimation of power spec - 
tral densities from periodograms (|Hauser and Peebleslll974l ). 
Basically, they estimate the Ci through the inverse Har- 
monical transform of a masked map that is then decon- 
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volved with geometrical kernels and corrected with a noise 
bias term. These techniques give unbiased estimates, and 
it has been shown they work successfully when applied to 
real data at high multipoles |jones et aL l2006l : IWu et aTl 
I2OO7I : IKuo et ali |2003 : lOunklev et al-lbood) . However, it is 
well known that at low multipoles they are not optimal since 
they provide power spectra estimates with error bars larger 
than the minimum variance. Several strategies for measur- 
ing d at low resolution have been developed and applied to 
CMB data with excellent results. These methods i nclude the 
Quadratic Maximum Likelihood (QML) estim ator l|Tegmarkl 
1 19971 : iTegmark and de Oliveira-Costal |2001^■ and di f ferent 



sampling technique s such as Gibbs (Jewell et al] 12004 
IWandelt et al. 2004 : Eriksen et all 12004 '). ada ptive impor- 
" 20091 7 and Hamiltonian (jTavlor et al.l 



tance (jBenabed et al 
I2OO7I ). 

In this paper 



we apply a 
mentation of the QML es ti mator 
iTegmark and de Oliveira-Costal I2OOII) 



parallel i mple- 



(|Tegmarkl 1 19971 : 

._ .. on WMAP 5 

year data. It is shown in l|Efstathioull2004l ') that at low mul- 
tipoles, the intrinsic variance introduced by this technique 
is lower than the variance introduced by the the pseudo-C^ 
methods. We also compute the conditional likelihood slices 
of the Ce multipoles and compare the results with those 
obtained by the QML. 

The largest angular scales of the CMB anisotropics 
probe the Physics of the early Universe. Therefore estimat- 
ing its power spectra with optimal methods is crucial. The 
temperature and polarization low multipole pattern contains 
information on the reionization process, dark energy and 
the geometry of the universe. Moreover, determining at best 
cosmological parameters as the optical depth r, has an im- 
portant impact on the amplitude and the spectrum of pri- 
mordial perturbations, and therefore on inflationary models. 
Optimal methods at low resolution are also required by the 
observed value of the temperature quadrupole, anomalously 
low within the ACDM dominant model which remains how- 
ever the simplest cosmological model preferred by data. 

In this paper we perform an analysis of the power spec- 
tra from low resolution maps of the five year WMAP data 
providing and errors bars for all the 6 APS from £ — 2 
to 32 computed by our QML and we show the robustness of 
our results with respect to iterative estimates. We also give 
the conditional likelihood slices computed with a pixel based 
likelihood code from £ — 2 to 10. The paper is organized as 
follows: in Section [5] we describe the main equation of the 
QML and in Section |3] we provide the basic expression for 
the likelihood. Section |4] is dedicated to the technical im- 
plementation with a brief description of the computational 
requirements both for the QML and for the likelihood codes. 
In Section [5] we report the results obtained from WMAP 5 
year data and finally in Section |S] we draw our conclusions. 
As a convention, all the objects that are defined in pixel 
space are represented with bold face symbols. 



2 THE QML METHOD 

The QML method for Power Spectrum Estima- 
tion (PSE) of CMB anisotropics was introduced in 
TegmarkI Il997l) and was ext e nded to polarization in 



in temperature and polarization x=(T,Q,U), the 
QML provides estimates Cf - with X being one of 
TT, EE, TE, BB, TB,EB - of the APS as: 
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with C = S{Ci) + N being the global covariance matrix 
(signal plus noise contribution, possibly extended to include 
residuals from foreground subtractions) and Ci is called fidu- 
cial power spectrum. It can be shown that the QML method 
is unbiased in finding the APS Cf'"^" of the map x: 



X,obs 



(4) 



It is also optimal, since it can provide the smallest error bars 
allowed by the Fisher-Cramer-Rao inequality. 



(5) 



where the averages are meant to be over an ensemble of 
realizations. 

We have verified by Montecarlo on simulated data that 
our implementation of the QML method (called BolPol, see 
Section U) leads to unbiased minimum variance results as 
from Eqs. (|4l[5]). We have also verified that Eq. (|4]) holds with 
Ci 7^ Ci''" , but the errors are no more minimum variance. 



3 CONDITIONAL LIKELIHOOD (OR 
PROBABILITY) SLICES 

Assuming that CMB anisotropics are Gaussian distributed, 
the likelihood for C MB temperature and p olarization power 



the likeuhood tor u ivm temperature and p ol; 
spectra is given bv lBond. Jaffe and Knoxl (jli 



1998)): 



exp 



ix*C- 
2 



/y/{2n)"det{C) 



(6) 



Tegmark and de Oliveira-Costal I2OOII ). Given a map 



where n is the dimension of the vector x (and depends on the 
number of pixels that are observed). It gives the probability 
to have the map x given the model (i.e. Ce) needed to build 
the signal covariance matrix. The conditional likelihood slice 
of Cf for given £ , X are obtained sampling Eq. (|6| at dif- 
ferent values of Cf keeping Cf fixed for £ ^ £ , X X . Our 
implementation in parallel fortran 90 of eq. (|6} is called Bo- 
Like (see Section|3]|. It is important to stress that the covari- 
ance matrix in x cither in BolPol or BoLike is treated with- 
out an y approx imation and not separated in T and (Q,U) 
as m (|Pageet al.ll2007l ). BoLike can also be employed as 
part of a likelihood code, to compute directly the posterior 
probability of a cosmological model at low £ directly in pixel 
space without relying on any approximation for the proba- 
bility distributions of the Cf . 

Other methods have been developed to compute the 
CMB likelihood for low resolution dataset, which can also 
be used to provide power spectrum esti mates. One such ap- 
proach is based on the Gibbs algorithm (jWandelt et al ]|2004l : 
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IChu et aljlioosh . which provides a computationally feasi- 
ble framework to sample the posterior probability of the 
Ce given the data, P(Cf|x). Specifically, this is achieved by 
repeatedly drawing samples from the (easier to compute) 
conditional distributions P(s|C^,x) and P{Ce,s), s being a 
signal only CMB map. A Blackwell-Rao app roximation can 
then be used to compute the CMB likelihood l|Rudiord et alj 
I2OO9I '). This me thod has been appli ed to the WMAP5 tem- 
perature data (|Dunklev et al 1 120091) and the WMAP 3 tem- 
perature plus polarization data l Eriksen et al.l [20071 ). One 
feature of the Gibbs approach is that it can internally per- 
form a parametric based component separation, and prop- 
agate model uncertain tes to power spectra and likelihoods 
jPickinson et al.ll200'i1 ). 

Another proposed approach to estimate the APS and 
likelihood at large angular scales is based on adaptive im- 
portance sampling. The idea is to model each single Ce dis- 
tribution with an inverse gamma distribution, finding appro- 
priate parameters. These distributions are then used to pro- 
vide an approximation to the global posterior. The method 
so far has been demo nstrated to work for intensity data and 
applied to WMAP5 (|Benabed et al..i2009i ). 



4 COMPUTATIONAL REQUIREMENTS 

BolPol is a fully parallel implementation of the QML method 
described in Section [2j written in F90. Since the method 
works in pixel space the computational cost increases as one 
considers smaller angular resolution for a given sky area, i.e. 
more pixels. This is the reason why we parallelize the QML 
code. The inversion of the covariance matrix C scales as 
O(A'^ix). The number of operations is roughly driven, once 
the inversion of the total covariance matrix is done, by the 
matrix-matrix multiplications to build the operators 'Eif in 
Eq. (|3]) and by calculating the Fisher matrix F^x' given in 
eq. @. The number of operations that are needed to build 
these matrices scales as 0{N'^i^^Npi^) and 0{NsideNpi^) re- 
spectively. The RAM required is of the order 0{A£ Np^^) 
where Ae is the number of C"^ (dC/dC^) (for every X) 
that are built and kept in memory during the execution 
time. 

Given these kind of scalings, it is clear that it is unreal- 
istic to run the QML estimator for all-sky maps of r esolution 
large r than Nside ~ 8 (in Healpix languag^J ( Go rski et al.l 
I2OO5I )) on a single CPU. To reach higher resolution we use 
the ScaLapack package and the BLACS routines. In this way 
it is possible to run BolPol on the WMAP data set with the 
resolution of N^ide = 16 on a supercomputer like BCX (at 
CINECA, Bologna, processor type: Opteron Dual Core 2.6 
GHz, with 4 GB per processor) in ~ 40 minutes using 64 
processors. 

The computational resources required by the likelihood 
evaluation in BoLike are driven by the calculation and the 
inversion of total covariance matrix. BoLike shares with 



^ http://healpix.jpl.nasa.gov/. For the reader not familiar with 
the Healpix notation, N^ide is related to the number of pixels 



A^pix by N, 



szde 



^ Note that N^ide = 16 is not the highest resolution that BolPol 
is able to consider. Currently BolPol is able to process maps of 



N. 



side 



: 32. 



BolPol the routines for the computation and inversion of 
the covariance matrix. 



5 DESCRIPTION OF THE WMAP 5 YEARS 
DATASET AND RESULTS 

In this Section we describe the data set that we have con- 
sidered. We use the ILC map in T, the foreground cleaned 
maps in {Q,U), the noise covariance matrix in (Q,U) and 
the masks at Nside = 16 publicly available at the LAMBDA 
web site 0. The temperature map is the ILC map at 
Nside ~ 16 smoothed at 9.8 degrees. We have added a ran- 
do m noise realiza t ion w ith variance of lfj,K^ as suggested 
in iDunklev et al.l l|2009t ) . Consistently the noise covariance 
matrix for TT is taken to be diagonal with variance equal 
to IfiK'^. We have explicitly checked that the estimates do 
not depend on the noise realization we added on the tem- 
perature map. The temperature map has been masked with 
a mask covering ^ 16 % of the sky and the monopole and 
the dipole have been subtracted from the observe d sky by 
mean s of the Healpix routine "remove_dipole" (Gors ki et al.l 
that works in pixel space. The polarization maps for 
Q and U are provided by the WMAP team at the same 
resolution Nside ~ 16. The inverse of the masked noise co- 
variance matrix for the polarization part is available at the 
same resolution. We have followed the procedure explained 
at the LAMBDA web site (see footnote 4) to obtain the 
direct noise covariance for the observed pixels. The used po- 
larization mask is larger than the temperature mask and 
covers ~ 26 % of the sky. The noise covariance matrices for 
TQ and TU are not provided and we set them to zero. In 
this work we use the WMAP 5 yrs publicly available data 
products to derive our spectral estimates. The uncertainties 
due to foreground cleaning cannot be taken into account 
since are not provided explicitly by the WMAP team. The 
QML estimates have been obtained by constructing the sig- 
nal covariance matrix by using fiducial Ces up to SA'^sidc and 
we self-consistently compute Cf up to SA'side- We show the 
QML estimates up to 2A'aide, which is a conservative choice 
for £-range, not to incur into discretization errors. 

We present the results obtained by our implementation 
of the QML estimator on the low resolution WMAP5 maps 
described in the previous section and those obtained by the 
WMAP team in Figs. 1-4, respectively for TT, EE, TE and 
BB spectra. In the top panels we show the estimates of 
BolPol with error bars (dark blue or red, see below) and the 
(pseudo-C^) estimates obtained by the WMAP team with 
error bars (light blue). The BolPol estimates in dark blue 
are obtained by u sing as fiducial spec trum the theoretical 
WMAP5 best-fit (|Dunklev et al.ll2009l ). a rACDM cosmo- 
logical model with n^h^ = 0.0227, n^h^ = 0.108, Ho = 
72 Akms'^ Mpc~^ , t = 0.089, = 0.961, As = 2.41 x 10~^ 
(at k — 0.002Mpc~^). Error bars loose dependence on the 
fiducial model by iterating the QML: we use the first run 
of BolPol^ to obtain another set of estimates with relative 



http://lambda.gsfc.nasa.gov/ 
* As fiducial spectra for the iterated BolPol run we use the first 
TT and TE BolPol estimates and leave the EE as given in the 
previous fiducial model when possible. The fiducial spectra of BB, 
TB, EB for the iterated BolPol run are set to zero. 
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Figure 1. Estimates of TT angular power spectrum from WMAP 5 year data at low resolution. Upper panel: BolPol estimates (dark 
blue diamonds) with error bars (dark blue), iterated BolPol estimates (red diamonds) with error bars (red), WMAP pseudo-C^ estimates 
(light blue cross) with error bars (light blue). Lower panel: differences between the sets of estimates in unit of sigma (same conventions 
as upper panel for the colors). 



error bar, which is plotted in red. The iterated estimates 
are always very close to those obtained with the WMAP5 
best-fit as fiducial model: this means that the QML esti- 
mated are sufficiently stable with respect to iteration. The 
same does not happen to the error bars: in particular the 
error on the TT quadrupole is substantially smaller, and a 
decrease of the error also occurs for ^ = 3 , 4 in temperature. 
Our estimate for the octupole in BB is consistent with zero 
and very different from the one obtained by WMAP (note 
however that the WMAP likelihood slice for CfJ^ is not 
anom alous as the WMAP pseudo-C^ estimate (|Nolta et al.l 

[iooi)). 

In the lower panels of Figs. 1-4 we show the dif- 
ferences between the sets of estimates in unit of sigma 
(same conventions as upper panels) . More precisely we show 
(^(jBoiPoi _ c^MAPy^ ^j^g BolPol estimates and a 

are given by: BolPol estimates with BolPol error bars for the 
no- iteration case (solid dark blue lines); BolPol estimates in 
the no-iteration case and WMAP error bars (dashed fight 
blue fine); BolPol estimates with BolPol error bars for the 
iteration case (solid red line). These lower plots show and 
quantify the differences among the considered sets of esti- 
mates. Our estimate of the TT spectrum is within la from 
the one computed by the WMAP team. For the other spectra 
we find larger differences (even 2a). The current joint anal- 
ysis is very different from the one performed by the WMAP 
team and it is obvious to expect larger differences for the 
more delicate spectra (where the Signal to Noise ratio is 
lower) . 



We list now the reduced 



X 



values for the iterated 



BolPol estimates from ^ = 2 to ^ = 32 with the fiducial 
input model: Xtt ~ 4.423 but excluding the quadrupole 
this value decreases to 1.079 recovering so the anomaly of 
the low quadrupole value of TT; for the other reduced 
values we find Xte = 0.785, Xee = 1-422, Xbb = 1-607. 

In Fig. 5, we plot the first and iterated BoiPo/ estimates 
and relative error bars (in blue and red, respectively) for 
TB and EB. We do not plot the estimates by the WMAP 
team since these are not provided in the LAMBDA site. Also 
for these parity-odd correlators the QML estimates are very 
stable with respect to the iteration. Note how the error bars 
in TB change (due to substantially different fiducial model 
in TT in the iterated run) , whereas those in EB do not (since 
the fiducial EE and BB spectra are mainly unchanged during 
the iteration). The TB nuU reduced for ^ = 2 - 23 is 1.34 
to be compared with 0.97 quoted for £ = 24 — 450 by the 
WMAP team. The EB nuU reduced for f = 2 - 23 is 1.14. 

At low multipoles, symmetric error bars, such as those 
provided by the Fisher matrix, are just an approximation, 
since we know that the likelihood for Ci is far from a sym- 
metric Gaussian. We therefore evaluate the conditional like- 
lihood slices for the six spectra from ^ = 2 to 10; we present 
these results in Figs. 6-7. As for the QML, we compute 
the slices on the WMAP5 best-fit (blue points) and on the 
BolPol estimates (orange points obtained with the same 
fiducial used for the iterated QML run). It is important to 
note that the peaks of the likelihood slices are very different 
for the two sets of conditionings: the quadrupole in tem- 
perature is the most striking example. We do not observe 
such dependence on the fiducial model in the estimates in 
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Figure 2. Estimates of EE angular power spectrum from WMAP 5 year data at low resolution. Upper panel: BolPol estimates (dark 
blue diamonds) with error bars (dark blue), iterated BolPol estimates (red diamonds) with error bars (red), WMAP pseudo-C^ estimates 
(light blue cross) with error bars (light blue). Lower panel; differences between the sets of estimates in unit of sigma (same conventions 
as upper panel for the colors). 



the QML method. As for the errors of the QML method, 
also the shape of the conditional slices depend on the fidu- 
cial model. The trend in the iteration is the same of the one 
observed for the QML: confidence levels shrink by insert- 
ing the iterated QML estimate as conditioning. In Figs. 6-7 
we also plot the QML estimate with error bars as reference 
(blue plus for the first BolPol run and orange for the it- 
erated one). The peaks of the likelihood slices are always 
consistent with the QML estimates within the error bars for 
both conditioning: however, an excellent agreement emerges 
when the QML estimates are used as conditioning. 

Few results deserve to be commented for their cos- 
mological importance. Our estimate (position of the peak) 
for the l{e + l)Cf^/(27r)|f=2 is 165 ^AT^ for the iterated 
QML (pixel base likelihood code). Basing on our pixel like- 
lihood code, the conditional likelihood (normalized to 1) for 
£{£+ l)Cf^/(27r)|fc2 is larger than 0.05 between 50 fiK^ 
and 1305 fiK^ . This las t value has to be com pared with the 
one given m Fig. 1 of iDunklev et all (|2009l ), in which we 
read that £{£ + l)Cf^ / {2n)\e=2 is approximatively less than 
1600 fiK'^ at the same confidence level. 

We also obtain constraints on BB. In Fig. [S]we show 
the conditional likelihood slice for the band powers, £{£ + 
l)Cf^/{2-K) binning from ^ = 2 to 6. Dark blue diamonds 
of Fig. |5]represent the likelihood slice we obtain sampling on 
band powers. We find £[£ + 1) Cf ^/(27r)|2-6 < 0.18/iif^ at 
95% of C.L. (this bound decreases to Q.lZfiK'^ at 95% of C.L. 
averaging over £ = 2 — 10). Light blue crosses in Fig. [S] rep- 
resent the conditional likelihood slice we obtain sampling on 
Ci. In this case we find £(£ + 1) Cf^ / {2-K)\2-.e, < Q.IQ^.K'^ 



at 95% of C.L. (this bound decreases to 0.11 ^,K^ at 95% 
of C.L. averaging over £ = 2 — 10). Some comments are 
in order: whereas it is expected to obtain different results 
sampling on different quantities, one may tentatively inter- 
pret the sampling on band powers more similar to a nearly 
scale-invariant spectrum of gravitational waves and the one 
on Cf^ similar to residual constant noise bias or a system- 
atic. Whereas we do not find evidence for a BB signal one 
may wonder about the different shapes of the conditional 
likelihood slices shown in Fig. |S] as a matter of fact, the 
quadrupole is the single multipole which contributes more 
to this discrepancy (see also the conditional likelihood slice 
for in Fig. 7. Excluding the quadrupole from the analy- 
sis, the conditional likelihood slice sampled on band powers 
yields an upper limit of £{£ + 1) Cf^ /{2%)\3-a < O.U^iK^ 
at 95% CL. While discrepancy is not statistically significant, 
the fact that it is mainly due to a single multipole might be 
pouting at a residual contamination in the WMAP 5 years 
map mostly affecting the quadrupole. 



6 DISCUSSIONS AND CONCLUSIONS 

We have performed a new estimate of the CMB angular 
power spectra at low multipoles from low resolution maps 
of the five years WMAP data. 

The QML estimates are found in agreement with the 
pseudo-Cf WMAP ones: the best agreement is for the TT 
spectrum and differences at the level of 2 — 3 a are found in 
EE and BB. 
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Figure 3. Estimates of TE angular power spectrum from WMAP 5 year data at low resolution. Upper panel: BolPol estimates (dark 
blue diamonds) with error bars (dark blue), iterated BolPol estimates (red diamonds) with error bars (red), WMAP pseudo-C^ estimates 
(light blue cross) with error bars (light blue). Lower panel: differences between the sets of estimates in unit of sigma (same conventions 
as upper panel for the colors). 
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Figure 8. Upper bounds on BB. y axis: Likelihood slice for BB. 
X axis i{e + l)Cf^/{2n) binning from £ = 2 to 6. Dark blue di- 
amonds represent the likelihood slice obtained sampling on band 
powers. Light blue crosses represent the likelihood slice obtained 
sampling on power spectrum. See also the text. 



origin of the diflFerence in the constraint in for £ = 2 — 6 
is partially due to this difference. 
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Whereas the QML is not used by the WMAP team, 
a pixel likelihood code is employed and results are given in 
iDunklev et al] (|2009l ). lNolta et all l|2009l ). The difference be- 
tween the conditional likelihood slices given here and those 
obtained by the WMAP team are due to our exact treat- 
ment of the covariance matrix: the approximation used by 
the WMAP team requires less computational resources but 
has the drawback of allowing negative values for conditional 
slices of EE and BB whereas ours are always positive. The 
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Figure 4. Estimates of TT angular power spectrum from WMAP 5 year data at low resolution. Upper panel: BolPol estimates (dark 
blue diamonds) with error bars (dark blue), iterated BolPol estimates (red diamonds) with error bars (red), WMAP pseudo-C^ estimates 
(light blue cross) with error bars (light blue). Lower panel: differences between the sets of estimates in unit of sigma (same conventions 
as upper panel for the colors). 
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Figure 5. BolPol estimates of TB (upper panel) and EB (lower panel) angular power spectra from WMAP 5 year data at low resolution. 
Dark blue symbols are for the not iterated case and red for the iterated case. 
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Figure 6. Likelihood Slices for TT, TE and EE from £ = 2 to 10 for the WMAP 5 year data at low resolution (i.e. nside = 16). Blue 
slices are for the not iterated case and the red ones for the iterated case. The blue plus represent the not iterated BolPol estimate with 
error bars (blue horizontal line) and the red plus the iterated ones. 
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